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We study the production of photons from a quark gluon plasma in local thermal equilibrium 
by introducing a non-perturbative formulation of the real time evolution of the density matrix. 
The main ingredient is the real time effective action for the electromagnetic field to 0{aem) and 
to all orders in as- The real time evolution is completely determined by the solution of a classical 
stochastic non-local Langevin equation which provides a Dyson-like resummation of the perturbative 
expansion. The Langevin equation is solved in closed form by Laplace transform in terms of the 
thermal photon polarization. A quantum kinetic description emerges directly from this formulation. 
We find that photons with k < 200 Mev thermalize as plasmon quasiparticles in the plasma on time 
scales t ~ 10 — 20 fm/c which is of the order of the lifetime of the QGP expected at RHIC and 
LHC. We then obtain the direct photon yield to lowest order in Oem and to leading logarithmic 
order in as in a uniform expansion valid at all time. The yield during a QGP lifetime t ~ 10 fm/c 
is systematically larger than that obtained with the equilibrium formulation and the spectrum 
features a distinct flattening for k > 2.5 Gev. We discuss the window of reliability of our results, 
the theoretical uncertainties in any treatment of photon emission from a QGP in LTE and the 
shortcomings of the customary S-matrix approach. 
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I. INTRODUCTION 



The quark gluon plasma (QGP) is a novel state of matter conjectured to be formed during ultrarelativistic heavy 
ion collisions and to have existed when the Universe was about 10 /isecs old. Ultrarelativistic heavy ion experiments at 
SPS-CERN, AGS-BNL, RHIC-BNL and the forthcoming LHC at CERN seek to create this state in collisions of heavy 
ions such as Pb and Au up to ^/s ~ 200 Gev/nucleon at RHIC and 5 Tev/nucleon at LHC . Current theoretical ideas 
suggest that a QGP in local thermodynamic equilibrium (LTE) is formed on a time scale of the order of ~ 1 fm/c after 
the collision when the partons in the colliding nuclei are liberated. Parton-parton collisions during a pre-equilibrium 
stage is then conjectured to lead to a state of local thermodynamic equilibrium that expands hydrodynamically and 
eventually undergoes a hadronization phase transition at a temperature of order ~ 160 Mev[l|. The experimental 
confirmation of a quark gluon plasma hinges on identifying observables that are directly associated with the properties 
of the QGP. Electromagnetic probes, namely photons and lepton pairs are considered to be "clean" since they only 
interact electromagnetically and their mean free paths are expected to be much larger than the size of the QGP. These 
probes are expected to leave the hot and dense region without further scattering, hence carrying direct information 
of the QGP0, 0, 3- These expectations led to an effort to obtain an assessment of direct photon emission and 
production of dilepton pairs from a thermalized QGP0-E3- Preliminary assessments concluded that direct photon 
emission from a QGP in LTE can be larger than electromagnetic emission from a hadronic gas|^||5|. 

The first observation of direct photon production in ultrarelativistic heavy ion collisions has been reported by the 
WA98 collaboration at SPS-CERN in ^ospb +2"^ Pb at = 158 Gev/nucleonQ. The WA98 data reveals an 
excess of direct photons above the expected background from hadronic decays in the range of transverse momentum 
Pt > 1-5 Gev in the most central collisions. These observations confirm the feasibility of direct photons as experimental 
probes for studying the formation and evolution of a QGP. 

A variety of fits of the experimental data with various theoretical models has been reportedly, however the results 
are inconclusive: models with or without QGP emission seem to fit the data in a manner qualitatively not very 
different from the fits based solely on hadronic 'cocktails' [lo|. 
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Photon data from ultrarelativistic heavy ion collisions at RHIC are forthcoming and we believe it is imperative to 
re-assess the current theoretical understanding of direct photon emission from a QGP in LTE for a reliable theoretical 
prediction, upon which may hinge the identification of a QGP. 

The most widely used approach to study photon production from a QGP in LTE is based on the S-matrix calculation 
of an inclusive transition probability per unit space-time volume with initial asymptotic states of quarks and giuons 
averaged with a thermal distribution. Recent studies in real-time have provided a detailed criticism of this approach 
on the basis that the S-matrix calculation explicitly assumes an infinite lifetime for a QGP in LTE and treats the 
uncofined quarks and giuons in the plasma as initial asymptotic states 15, 16J. The underlying assumptions in the 
widely used S-matrix approach, namely, treating the medium as of infinite lifetime and spatial extent are manifestly 
inconsistent with the experimental and physical situation in which the state prior to the collision is described by 
hadrons and the ensuing QGP, if formed, only lasts during a (proper) time of order lOfm/c and its spatial scale is 
~ 10/m. There are preliminary studies of the influence of finite size effects on photon emission from the plasma [T^ 1 18| | 
and more recently for dilepton emission from the hadronic gasj^^. Studies of the finite lifetime effects on the photon 
production yield were reported in references ITfil|. The results of these studies point out the importance of non- 
equilibrium real time processes whose contribution is subleading in the limit of infinite lifetime but that during the 
finite lifetime yield contributions of the same order of or larger than those found in the S-matrix approach. 

Studying particle emission from a QGP is fundamentally different from a scattering experiment. In a scattering 
experiment the beam consists of physical particles, namely asymptotic "in" states, while in a QGP the quarks and 
giuons exist as a deconfined state of matter as a transient state. The "in" and "out" states are hadrons not quarks. 
The QGP thus emerges as an intermediate, transient state, and treating it as stationary source of electromagnetic 
radiation as is explicit in the S-matrix approach should be taken, at best, as an approximation. 

The fundamentally correct calculation of photon production in an ultrarelativistic heavy ion collision must neces- 
sarily begin from an in state of nuclei formed by confined quark and giuons. These bound states evolve from the 
infinite past with the full Hamiltonian continuing the evolution through the collision and deconfinement process, 
through the thermalization and formation of a thermal QGP, through hadronization and eventually through final 
state interactions after freeze out. The photons produced during all of these processes will be measured as out states 
in the detector. Photons are produced at every stage, and the total number of photons can be calculated from the 
S-matrix if all of the processes are included in the dynamics with the inclusive out states being hadrons. Of course 
this is so far an unsolved theoretical task, and the problem emerges when one attempts to isolate individual stages 
of the dynamics and treat these stages, each of which has a finite time duration, in terms of an S-matrix calculation. 
For example, current estimates neglect the photons produced prior to the collision, treats the pre-equilibrium stage 
via parton cascade models, and ignores the photons from this stage in the calculation of emission from a QGP, which 
itself is treated as a steady state 15j. 

The potential identification of the formation and evolution of a QGP from electromagnetic probes from forthcoming 
RHIC data requires a reliable theoretical understanding of the phenomena. 

In particular a recent detailed study of photon production in real and finite time from a QGP in LTE reveals 
that the spectrum of photons emitted during the finite lifetime of the QGP depends on the initial state at the 
time of thermalization and of the pre-equilibrium stage. This dependence is more pronounced for large momenta 
/c > 4 — 5 Gev. Such conclusion is in agreement with those obtained in ref.|l3| wherein a sensitivity of the large kx 
part of the spectrum to initial conditions was also found. The real time studies of photon production from a QGP in 
LTE with a finite lifetime in rcfs. 15, 16J all reveal a flattening of the spectrum at large momentum, the precise value 
of the momentum at which the spectrum flattens being sensitive to the pre-equilibrium stage which determines the 
initial state of the plasma at the thermalization hypersurface. 

Perhaps coincidentally the WA98 data 14] displays a flattening of the spectrum for pt ^ 1-5 Gev. 

Windovir of opportunity: Before a calculation of the direct photon yield from a QGP in LTE is attempted, 
it is important to establish the regime of experimental relevance and of theoretical reliability of any prediction 
of the yield based on LTE. The region of soft photons k < 100 — 200 Mev is experimentally complicated by the 
enormous background of photons from neutral pion decay (produced profusely in ultrarelativistic heavy ion collisions), 
bremsstrahlung from final state interactions etc. Assuming that LTE is established and maintained by quark and gluon 
collisions and assuming a coUisional mean free path of order of 0.1 < A < 0.3 fm, photons with momenta fc > 3 — 5 Gev 
will likely probe scales shorter than the mean free path where the LTE approximation breaks down. Large departures 
from hydrodynamic behavior entailed by LTE is revealed in recent elliptic flow data on the parameter f9(pT)pol l2l| 
for pt > 2 Gev. While this data indicates the breakdown of (ideal) hydrodynamics for hadrons, a similar breakdown 
is expected on physical grounds for hard electromagnetic probes when the Compton wavelength of the probe is smaller 
than the scale of the mean free path. Thus the reliability of a calculation of hard particle emission from a QGP in 
LTE must be re-assessed. 

Furthermore the spectrum of high energy photons for /c > 4 — 5 Gev is sensitive to the initial conditions an the 
pre-equilibrium staeefTsl Il5|. Originally the pre-equilibrium stage was modelled as a parton cascade|i22j. however 
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more recently a different picture is emerging for the description of the pre-equihbrium stage based on color glass 
condensates . Clearly the current understanding of the pre-equilibrium stage is still a matter of ongoing study. 

Thus from the experimental point of view the relevant range of "clean" photons is probably for k > 100 — 200 Mev, 
and from a theoretical point of view, the reliability of any calculation of the yield based on the assumption of LTE 
is probably suspect for A: > 4 — 5 Gev with the added uncertainty of dependence on the pre-equilibrium stage for 
momenta larger than this range. Therefore, in this article we will focus our study to the range 200 Mev < A: < 5 Gev 
commenting on further uncertainties emerging from our study for the high energy region. 

Goals of this article: In this article we study non-perturbative aspects of the real time dynamics of photon 
production and propagation with the goal of obtaining a deeper understanding of the direct photon yield during the 
finite lifetime of a transient QGP in LTE. In particular we focus on establishing a real time description of photon 
production and obtaining a theoretical prediction of the yield and the spectrum by including processes that are 
not included in the usual S-matrix calculation even at lowest order in aem- These processes lead to subdominant 
contributions in the asymptotically long time limit, but their contribution during the finite lifetime of the QGP is of 
the same order of or larger than those extracted solely from an S-matrix analysis. 

We begin by obtaining the time evolution of an initial density matrix directly from the effective action for the elec- 
tromagnetic field exact to order and to all orders in the strong coupling a^. This formulation makes manifest the 
connection between photon production and the stochastic nature of photon emission and propagation in a thermal- 
ized plasma. The expression for the photon production yield obtained from this description reproduces the S-matrix 
results in a strict perturbative expansion in for a steady state QGP. Moreover, when the full evolution is taken 
into account this formulation includes the dynamics of photon propagation in the medium and of thermalization. 
Furthermore, this formulation reproduces the results of kinetic theory, highlighting the limitations of the equilibrium 
approach to photon production. 

Summary of main results: 

• We provide a formulation of direct photon production in real time by obtaining the effective action for the 
electromagnetic field up to lowest order in aem and in principle to all orders in as- Integrating out the quark 
and gluon fields to obtain the effective action leads to the description of the production and propagation of pho- 
tons in a thermal bath. The properties of the QGP thermal bath are determined by current-current correlation 
functions and a stochastic gaussian colored noise that obey a generalized fluctuation-dissipation relation com- 
pletely determined by the thermal photon polarization. The time evolution of the photon distribution function 
is determined by the solution of a classical stochastic and non-local Langevin equation. We explicitly solve the 
non-local Langevin equation in closed form by Laplace transform in terms of the thermal photon polarization. 
The resulting evolution represents a Dyson-like resummation of the naive perturbative expansion and provides 
a uniform expansion in valid at all times and which includes non-perturbative aspects. We obtain the 
photon yield and the spectrum during the lifetime of the QGP expected at RHIC (~ 10 fm/c) up to leading 
order in aem and to logarithmic order in a^. 

• An important aspect that emerges from our study is that intermediate energy photons with k < 400 — 500 Mev 
propagate in the QGP plasma as quasiparticles, and in particular photons with k < 200 Mev thermalize with the 
plasma on time scales of the order of '--^ 10 — 20 fm/c which is of the order of the lifetime of the QGP expected 
at RHIC. 

• We predict the photon spectrum for the range of momenta 200 Mev < < 5 Gev. The yield is calculated over a 
time scale compatible with the expected lifetime of a QGP at RHIC or LHC ^ 10 fm/c. In the intermediate region 
200 Mev ^ k <2 Gev the spectrum is similar to that obtained by the S-matrix formulation but systematically 
larger by a factor that ranges between 2 — 4. The spectrum flattens at an energy scale k ~ 2.5 Gev becoming 
dramatically larger than the yields obtained previously with the S-matrix formulation. 

The article is organized as follows: in section II we obtain the real time effective action to lowest order in aem 
and in principle to all orders in as- This effective action manifestly establishes contact with the stochastic nature of 
photon emission from a plasma in thermal equilibrium. In section III we obtain the photon distribution function in 
real time in terms of the solution of a Schwinger-Dyson equation which includes the photon self-energy to lowest order 
in aem and in principle to all orders in as ■ This solution provides a Dyson-like resummation of the naive perturbative 
expansion. In section IV we show that the S-matrix result emerges in the strict perturbative limit. In this section 
we also make contact with a kinetic description of photon production which manifestly includes the dynamics of 
thermalization. In section V and VI we study the non-perturbative aspects and obtain the main results of this article. 
Section VII summarizes our conclusions and presents further questions. 
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II. THE REAL-TIME EFFECTIVE ACTION 



In refs. jld l24| a manifestly gauge invariant formulation of the time evolution of an initial density matrix has been 
described. We will follow this treatment as it guarantees that the results are completely gauge invariant. The gauge 
invariant Hamiltonian is given bv [l5llT^l24| 



H = HQcom + j d?x- [El + B^) + e J d^x 3 ■ At + H.^ui , (H.l) 

where Hqcd [^] is the QCD Hamiltonian in absence of electromagnetism but in terms of the gauge invariant (under 
abelian gauge transformation) quark ficld(s) ^E" and the subscript T refers to transverse components. We have extended 
the fermion content to Nf flavors and the charge of each flavor species in units of the electron charge is included in 
the current, namely 



Nf 

3 = J2-^^l^^, (n.2) 

The instantaneous Coulomb interaction can be traded for a gauge invariant Lagrange multiplier field which we call 
(which, however should not be confused with the time component of the gauge field), leading to the following 
Lagrangian density 



1 ^' 

^ = ^QCD + /:o,em - J°A" + J • At ; /:o.em - [(^^At)' + (VA0)2] . ^ ^ fi ^ (H.S) 

We will study the time evolution of the number of physical transverse photons and the expectation value of the 
(gauge invariant) transverse gauge field in the linearized approximation as an initial value problem. 

In ref.^^ correlated initial states in which quark states are dressed by the electromagnetic field were studied in 
detail. One of the important conclusions of that study is that such an initial density matrix cannot correspond 
to a QGP in LTE under the strong interactions. This is a consequence of the fact that correlated quark-photon 
states necessarily require states constructed with the quark-current operator, which does not commute with Hqcd- 
Furthermore the results of this reference suggest that the details of the initial pre-equilibrium stage become manifest 
in the spectrum for fc ^ 4 Gev. The dependence of the high energy part of the spectrum on initial conditions has also 
been studied in^^ with similar conclusions. 

Here we study the simpler case of an uncorrelated initial density matrix to highlight the non-perturbative aspects, 
with the understanding that a precise description of the high energy region will necessarily require a firmer knowledge 
of the initial state and the pre-equilibrium stage. 

Hence we propose the initial density matrix to be of the form 

p{ti) = pQCD-L ® PAT,i (II. 4) 

The initial density matrix pQCD,i describes a quark-gluon plasma in (local) thermodynamic equilibrium at a tem- 
perature T — namely 



PQCD,,: = e-^^«^- (H.5) 

And ,i is diagonal in the Fock representation of free field quanta of the transverse gauge field with an initial 
distribution of quanta. It is important to emphasize that any other initial density matrix that mixes photons and 
quarks will not commute with iJQcu'lS], hence cannot describe a QGP in equilibrium. 

The vacuum state is represented by PAT,i — with |0) being the vacuum state of the field At{x). 

In the field basis (Schroedinger representation) the matrix elements of pAT,i ^re given by 

(AtIpatA^t) ^ PAtA^t;At) (II-6) 
The time evolution of the initial density matrix is given by 



-iH{tf-U) 



piu] 



^iH{tf-ti) 



(H.7) 
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where H is the total Hamiltonian given by eq. This particular choice of initial state will introduce transient 

evolution, however the long time behavior should be insensitive to this initial transient. 

Furthermore, we point out that it is important to study the initial transient stage for the following reason. Photons 
are produced in the thermal bath since the initial density matrix does not commute with the total Hamiltonian 
(non-equilibrium) and they propagate in the medium as "quasiparticles" . By studying the initial transient dynamics 
after photons and quarks are coupled we can address the question of the dynamics of the formation and propagation 
of the quasiparticle which will be studied in detail in section ljVI|l . 

Our goal is to obtain the real time evolution of the number of photons produced as well as that of the equation of 
motion for the expectation value of the (gauge invariant) transverse gauge field to lowest order in aem but in principle 
to all orders in ■ The equation of motion for the expectation value of the transverse gauge field will yield real time 
information on the formation and propagation of transverse plasmon quasiparticles {see section IVll . 

The strategy is to obtain the real-time effective action for the transverse gauge fields by integrating out the quark 
and gluon degrees of freedom to lowest order in aem but to all orders in as- In this manner, the QGP in LTE is 
treated effectively as a thermal bath. 

The calculation of correlation functions is facilitated by introducing currents coupled to the different fields. Fur- 
thermore since each time evolution operator in eqn. (|II.7|I will be represented as a path integral, we introduce different 
sources for the various fields for forward and backward time evolution operators, referred generically as respec- 
tively. The forward and backward time evolution operators in presence of sources are U{tf,ti]J'^), U^^{tf,ti,J~) 
respectively. 

In order to avoid cluttering of notation let us collectively denote by x the quark and gluon fields which will be 
integrated out 

In what follows we will ignore the Lagrange multiplier field Aq since the Couloumb interaction will be irrelevant to 
leading order in aem- 

The non-equilibrium generating functional then is given by 

z[]+,n = T:YU{^M,J^)p(ti)u-\^,t,,j-) = 

J DAt,, J DAt^,paU^t.,;At,,) J I?A± J i?;^±e^^[^-'^*^-'-T^-'^l (II.8) 

where 

dtd^X [Co.em{A+) + J+A+ - /:o,e,„(Ay) - J^Ay] + 

J^d\[CQCD{x) + Jx + e3 -At] (II.9) 

and C describes a contour in the complex time plane as follows: from ti to +00 (forward) the fields and sources are 
AJ, x"*", (where x collectively represents both quarks and gluons), from -l-oo back to ti (backward) the fields and 
sources are A^, x~, and from ti to ti — i(3 (Euclidean or Matsubara) the fields and sources are At = 0; x'^i Jx- 
Along the Euclidean branch the interaction term vanishes since the initial density matrix for the quark and gluon 
fields, generically denoted as x is assumed to be that of a QGP in thermal equilibrium. One can in principle consider a 
correlated initial state of a QGP and electromagnetic fiuctuations but such state does not describe a QGP in thermal 
equilibrium and such initial density matrix would not commute with Hqc d , contrary to the usual assumption of an 
equilibrated QGPffSj. 

We seek to obtain the real time effective action for the transverse gauge field. Therefore we integrate out (trace 
over) the quark and gluon fields. After carrying out the trace over the QCD degrees of freedom, the remaining path 
integrals over At with the real time effective action have boundary conditions on the field At given by 



A+(f,t = t,) = AT,i(f) ; Ar^{x,t = t,) = A'j,^^{x) (11.10) 

The real time effective action for At is obtained by treating the quark and gluon fields as a bath performing the 
path integral over the quark and gluon degrees of freedom, namely, tracing over the bath degrees of freedom, leads to 
the influence functional ji25j| for A^. 

The initial density matrix for the At field will be specified later as part of the initial value problem. 

To lowest order in aem but to all orders in as the real time effective action for At is obtained as follows. Insofar 
as the path integral over the QCD degrees of freedom is concerned, At is simply a background field, hence 



eJ- A3 



)qcd 



(11.11) 
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expanding the expectation value in powers of e 



(11-12) 

where we have used {3)qcd = and the connected current-current correlation function {3^ 3^) qc d is in principle to 
all orders in as- 

We introduce the spatial Fourier transform of the quark current J as 



mi) 



1 



rf-*a;e"=-^J(f, t) 



with f2 the quantization volume, in terms of which we obtain following the correlation functions 

e^m\t)hn{-k-.t')) = e\iT{k-t)i+{-k-t')) = T>^{k;t-t') = T.i+{k-t,t') 

e^(3m{-k;t')M{kt)) = e^j+ik;t)j-{-k;t')) = T.<^{k-t-t') = T,t-{k-t,t') = 

e^{T ji{k;t)jmi-k;t')) = T>^{k;t - t')Q{t - t') + T.<^{k-,t ^ t')Q{t' ^ t) =^++ 

e\f ii{k;t)im{-k;t')) = T.>^{k;t- t')Q{t' -t) + T.Uk;t- t')Q{t - t') = ^i-{k;t,t') 



(11.13) 



(11.14) 
(11.15) 
(11.16) 
(11.17) 



where the expectation values are in the equilibrium thermal density matrix of the QGP, which results in translational 
and rotational invariant correlation functions that are only functions of the time differences and the superscripts ± 
refer to the forward (+) or backward (— ) time branches. The symbols T and T refer to time and antitime ordering 
respectively. These correlation functions are not independent, they obey 



S+„+(fc;t,0 + Er„-(fc;t,0 - ^i+{k;t,t') - ^t-{k;t,t') = 



(11.18) 



It is convenient at this stage to separate the transverse and longitudinal components of the polarization tensor 
^im{k; t — t') as follows 



^im{k;t-t') = Vimi'k)^T{k;t-t')+-kikm^L{k;t-t') 
^im(k) = Sim - k/k^ 



(11.19) 
(11.20) 



The non-equilibrium real time effective action in terms of the spatial Fourier transforms of the fields and correlation 
functions to lowest order in agm but to all orders in as is therefore given by 



iSgf AA+,A^] = Y^IU dt[A+^{t).A+.^{t)-k'A+^{t).A+.^{t) 



-Az (t)A-^ (t) + k^Az (t)A-^ (t) 



dt j dt' [A+^{t) ■ A+. ^(t')S++(fc; t-t')+ AZ^it) ■ ^(i')S^-(fc; t - t') 

~A+Jt) ■ A:.^(t')S+-(fc;t - t') - Ar^(t) . A+-^(t')S?+(^; t - t')] } (11.21) 

As it will become clear below, it is more convenient to introduce the following Wigner center of mass and relative 
variables for the transverse gauge field 



A{x, t) = ^ {At (^) t) + Arp {x, t)) ; a(x, t) = A;^ {x, t) — Aj, (x, t) 
and the Wigner transform of the initial density matrix for the transverse gauge field 



(11.22) 



(11.23) 
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The initial conditions on the At path integral given by ljII.10|l translate into the following initial conditions on the 
center of mass and relative variables 



A{x, i = 0) = A{x) ; a(f , t = 0) = a.. 



(11.24) 



The center of mass variable plays an important role: its expectation value is the mean field, since the expectation 
values of A^{x,t) coincide, 



{A+{x,t)) ^TrAT{x,t) p^Trp At(x,<) = {A:j,{x,t)) 



(11.25) 



In terms of the spatial Fourier transforms of the center of mass and relative variables (|11.22() introduced above, 
integrating by parts and accounting for the boundary conditions (|II.24p the non-equilibrium effective action llll.21|) 
becomes: 



iSeff[A,a] = I dtY^^-i d_j:{t) ■ (A^{t) + k^A^it)^ } 



dt I dt' { -d_j:{t) ■ %(t')/CT(fc; t~t')+ d^^{t) ■ Az{t') iY.^{k; t ~ t') 



+ / d'^xdiix) ■ A{x,t:^ 0) 



(11.26) 



where the last term arises after the integration by parts in time, using the boundary condition (|II.24|1 . There is no 
contribution from the t — > oo limit since the fields Aj^lt) at non-zero temperature will vanish at asymptotically long 
time. The kernels in the above effective Lagrangian are given by 



lCT{k;t-t') = ^[J:^{k;t-t') + ^:f{k;t-t')] 

iE^{k;t-t') = [T,^{k;t-t') - Y,^{k;t - t')] e(t-t') = iT,Tik;t ~ t')e{t - t') 



(11.27) 
(11.28) 



where we have used the fact that the polarization is a function of the modulus of the wavevector by rotational 
invariance. The photon polarization is computed to lowest order in and in principle to all orders in ag. 

The quadratic part in e^^''ff^-^''^^ in the relative variable a can be written in terms of a stochastic noise variable ^ as 



1 



exp 



dt / dt'd_^{t) ■ d^{t')JCT{k\ t - t') 



2?Cexp<i-i / dt I dt' Cj:/t)IC^Ht - t')^.kjM') 



+z J dt ^_fe/t) •«£,,(*) I (11.29) 

The non-equilibrium generating functional can now be written as 

Z = I DA^ I D£, I VAVdV^ W{A^■,£i)D ^^g^/ (u go) 



exp < — i j dt a_j:^{t) 



Aj:,{t) + eAj:^{t) + / dt' E#(fc; t - t')Aj:,{t') - ^^.(i) 



where the noise probability distribution function is given by 



PKl = cxp|-i jdtj dt'J2 



(11.31) 



The functional integral over Si can now be done, resulting in a functional delta function, that fixes the initial 
condition Arix^t = 0) = AT.iix). 

Finally the path integral over the relative variable a can be performed leading to a functional delta function and 
the final form of the generating functional is given by 



Z= DA^DS, WiA;£^)VAV^ ^K]%;(5 



Aj:i{t) + k^Aj:i{t)+ dt' ^T{k;t ~ t')A^ ^{t') ~ ^{t) 
Jo ' ' . 



(11.32) 
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with the initial conditions on the path integral on given by 

A{x, t^Q)= A^{x) ; A{x, t = 0) = £,{x) (11.33) 

and we have used the definition of (fc; t — t') in terms of ST(fc; t — t') given in equation l|II.28|l . 

The meaning of the above generating functional is the following: in order to obtain the correlation functions of 
the center of mass Wigner variable A we must first find the solution of the classical stochastic non-local Langevin 
equation of motion 

Aj:i{t) + k^Aj:i{t)+ I dt' ST(fc; t - = ,W 

Jo 

A{S, t = 0) A{x) ; A{x, t 0) = £^{x) (11.34) 



for arbitrary noise term ^ and then average the products of A[S,] over the stochastic noise with the probability 
distribution given by (CUni, and finally average over the initial configurations Ai{x) ; £i{x) weighted by the 
Wigner function W{Ai; £i), which plays the role of an initial semiclassical phase space distribution function. 



Calling the solution of pi.34|l A^ ^{t; Ai; Si) , the two point correlation function, for example, is given by 
{A_j:iit)Aj:^^{t')) = Z'' J V[^]V[^] J DA, J D£, W{A;£,) ^(i; A; ^1) ^^/t'; A; 



xS 



Aj:,it) + k'Aj:,{t)+ dt' ^T[k-t-t')A~^^,{t')-t^^^{t) 







(11.35) 



We note that in computing the averages and using the functional delta function to constrain the configurations of 
A to the solutions of the Langevin equation, there is the Jacobian of the operator 



5{t~t') 



'2 



2 



T.^{k]t--t') (11.36) 



which, however, is independent of the field and cancels between numerator and denominator in the averages. 

This formulation establishes the connection with a stochastic problem and is similar to the Martin-Siggia-Rose(2^ 
path integral formulation for stochastic phenomena. There are two different averages: 

• The average over the stochastic noise term, which up to this order is Gaussian. We denote the averages over 
the noise with the probability distribution function P[^] given by eqn. Ijll.Sip as 

Since the noise probability distribution function is Gaussian the only necessary correlation functions for the 
noise are given by 

= , ((%/i)efc,,,(i')» = ^y(k) ICrik^t- t')S'{k + k') (11.38) 
and the higher order correlation functions are obtained from Wick's theorem. 

• Average over the initial conditions with the Wigner distribution fmiction W{Ai;£i) which we denote as 



JDAJD£, W{A;£,) 

The initial Wigner distribution function requires precise knowledge of the pre-equilibrium stage. It allows to 
include initial correlations of quarks and photons, namely "entangled states" but as discussed in detail in ref.|l5j| 
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such initial density matrix will not commute with the Hqcd Hamiltonian and will not describe a state in LTE 
under the strong interactions. All of the theoretical uncertainties about the initial state prior to equilibration 
are encoded in the initial density matrix or alternatively in the initial Wigner distribution function yV(^i;£i). 
In ref.^^ a particular state with initial correlations reflecting a pre-equilibrium stage was modelled. While we 
can consider such modelled initial state, in this article we focus on extracting non-perturbative aspects in the 
simplest and cleanest scenario, that of an initially uncorrelated Gaussian state for photons. The numerical study 
performed in ref.O^ revealed that initial state preparation on time scales of 0(1 fm/c) modifies the photon 
spectrum for hard momenta k > 4 Gev. This is roughly the scale of momenta at which the assumption of 
LTE breaks down in any case because the photon probes scales of the order of the mean free path for quark- 
gluon collisions. Thus studying an uncorrelated initial state, which is compatible with all previous S-matrix 
calculations is phenomenologically relevant for photon momenta fc < 4 — 5 Gev. 

Therefore in what follows we will consider a Gaussian initial Wigner distribution function with the following 
averages: 



i.k i. — k 



i.k i, — k 



7"'''(k) 
2k 



[1 + 2A4-] +.4"- - 



V'^''{k)-[l + 2Af,]+£^^ £l^.^ 



A'^ _^ + _^A'^ ^ ^ A"^ ~ _ - + c.c 



(11.40) 

(11.41) 
(11.42) 



where A/fe is the initial distribution of photons. The indices a, b refer to vector components, while the label i 
refers to the initial conditions (III. 331). 



Thus averages in the time evolved full density matrix are therefore given by 



(•••)-((•••» • (11.43) 

The noise field emerged as an auxiliary variable to represent the quadratic contribution from the relative variable in 



eqn. Ijll.26(l . While the full microscopic dynamics is completely Hamiltonian and therefore deterministic, integrating 
out the QGP degrees of freedom leads to a reduced density matrix and the ensuing coarse grained dynamics for the 
photon field, just as in the microscopic treatment of Brownian motionpsj]. 

A. Relation between E and )C : Fluctuation and Dissipation 

From the expression (|11.28|l for the photon polarization which is determined by averages in the equilibrium density 
matrix of QCD we now obtain a dispersive representation for the kernels ICT{k;t — t') and S^(fc;t — t'). This is 
achieved by explicitly writing the expectation value in terms of energy eigenstates of the bath (quark and gluon fields) 
introducing the identity in this basis and using the time evolution of the Heisenberg field operators to obtain 

/oo 
du a>{k;uj) e''^^*-*"> (11.44) 
-OO 

/OO 
CT^(/c;w) e^"(*-*') (11.45) 
-OO 

with the spectral functions 

a>{k;u;) = e^^g^i^^ V e"^^" (n|j,(fc, 0)|™) (m|j, (-fc, 0)|n) 5(c. - (i?„ - i?™)) (11.46) 

2^QCD 



^T{k\^) = e2:^#^Ve-^^"(n|j,(-fc,0)|m)(m|j,(fc,0)|n)5(c.-(i?„,-i?„)) (11.47) 

^^QCD — 



m.n 
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where Zqcd is the QCD equiUbrium partition function, the states \n > are exact eigenstates of Hqcd and we have 
used rotational invariance. Upon rclabcUing m ^ n in the sum in eq. (|11.47(l we find the KMS relationpTl 



{k;u}) — CT^ (fc; — oj) = e'^'^cr^ (fc; uj) 



(II. 



where we have used parity and rotational invariance in the second line above to assume that the spectral functions 
only depend of the absolute value of the momentum. 

Using the spectral representation of the 8(t — t') we find the following representation for the retarded photon 
polarization given by eqn. ^(Hl^ 



with 



(11.49) 



duj 



uj — ko + ie 

Using the condition (|11.48|) the above spectral representation can be written in a more useful manner as 



(11.50) 



t^ik.ko) 



ImSP(fc,tj) 



duj- 

-oo uj- ko + ie 



ImS^(fc,w) = 7ro->(fc;w) [e'^" - l] 



(11.51) 
(11.52) 



clearly ImE^(fc; w > 0) > 0. Eq. (|II.48(I entails that the imaginary part of the retarded photon polarization is an odd 
function of frequency, namely 



ImSp(fc,w) = -ImS^(A:, -w) 

which is manifest in eq. pi.52|l . 

The relation (|11.52p leads to the following results which will prove to be useful later 



(11.53) 



a}^{k;uj) — — ImE|?(fc, w) n((jj) 

TT 

a^{k;uj) = -ImE^(fc;w) [1 + n(tj)] 

TT 

Similarly from the eqs. Ijll.27p and 111.4411 . (|II.45|I and the condition (|II.48|) we find 



(11.54) 
(11.55) 



J-oo 27r 

trik-ka) = 7ra>(fc;fco) [e^'^' + l] 



whereupon using the condition ljll.48|l leads to the generalized form of the fluctuation-dissipation relation 



(11.56) 
(11.57) 



iCT{k\ fco) = ImE^(fc; fco) coth 



/3fco 



(11.58) 



Thus we see that ImE^ (fc; fco) ; JCxik; fco) are odd and even functions of the frequency respectively. 

For further analysis below we will also need the following representation for Sr(fc; t — t') introduced in eq. HII.28P 



/oo -poo 
^^^{t-t•) _ ^< duj = _ e''^(*"*')lmE^(fc; Lj)duj (11.59) 

-OO ^ J — OO 
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whose Laplace transform is given by 



ET(fc;s)= / dte-'*XT{k;t) = -- r ^-^^^^^^^dcj (11.60) 

TT J_r^ UJ + IS 







This spectral representation, combined with 111.51(1 lead to the relation 

E^(/c;w) = ST(fc;s = + (11.61) 

III. THE PHOTON DISTRIBUTION FUNCTION IN REAL TIME 

The solution of the Langevin equation (|II.34p can be found by Laplace transform. Defining the Laplace transforms 



^kA^)^ dte--^s/0 ; ^ / rfte-%,,(t) (IILl) 

along with the Laplace transform of the photon polarization given by eqn. I(IL60() we find the solution 

£1 + sAl , + ij: lis) 
, ^ _M ''''^ '^-'^ (IIL2) 

where 



4^ = 7^ / d'xe*'~MKx) ; = / d'xe^-'^Elix) (IIL3) 



and the superscript i refers to the initial conditions (|II.33|) . The solution in real time can be written in a more compact 
manner as follows. Introduce the fundamental solution fk{t) of the equation of motion 



fk{t) + fk{t) + I dt' Erik; t - t')fk{t') = (III.4) 

obeying the initial conditions^ 



Its Laplace transform is given by 



/fc(t = 0)=0; /fe(t = 0) = l (IIL5) 



Ms) = , (nL6) 

+ fc^ + ijT(k; s) 



which is recognized as the Laplace transform of the full transverse photon propagator. 
The fundamental solution fk{t) is found by the inverse Laplace transform 



/,(t) ^ / ^ €L (IIL7) 

'c 2TTis^ + k^ + Y.Tik;s) 



^ The lower limit in the integral t = simply reflects the choice of the initial time. If an arbitrary initial time is chosen to, the lower 
limit becomes to, since the T,T{k;t — t') is manifestly time translational invariant, the solution of the equation of motion is a function 
of t — to- 
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where C stands for the Bromwich contour, parallel to the imaginary axis in the complex s plane to the right of all the 
singularities of /(s) and along the semicircle at infinity for Res < 0. The singularities of /(s) in the physical sheet 
are isolated single particle poles and multiparticle cuts along the imaginary axis. Thus the contour runs parallel to 
the imaginary axis with a small positive real part with s = iuj + e ; — cx3<cl;<oo and wraps around returning parallel 
to the imaginary axis with s — iuj — e;oo>uj> —oo, with 6 = 0^. From the spectral representations ljII.52III.60|) ) 
one finds that ET(fc, s = ioj zt e) = ReSj;(fc, oj) ± ImS.p(fc, ui) and using that ImE^(A:, us) = — ImS^(fc, —ui) we find 



fk{t) = 



dco 



sm{ujt) Iml]^(fc;a;) + 2u!e 



- fc2 - Re5]^(fc;w) + lm.t^{k;uj) + 2uje 



We have kept the infinitesimal 2uje ; e ^ 0+ to highlight the possibility of isolated quasiparticle poles in the case of 
vanishing imaginary part of the polarization. 

The initial condition /fe(i = 0) = 1 results in the following sum rule 



TT 



ImE^(fc;tj) + 2uje 



ReE^(fc,: 



ImE^(fc;w) + 2uje 



(III.9) 



In terms of the fundamental solution fk{t) given above, the solution of the Langevin equation (|II.34(I in real time 
is given by 



Aj: ^{t; A-£,- = A'f: ^ h{t) + f i , fk{t) + / .fk{t - t') ^i^i{t')dt' 



(III.IO) 



where the superscript i refers to the initial conditions ljll.33|) . 

This solution in real time represents the resummation of the Dyson series. 



A. Number operator: 



We define the photon number operator to be given by 

A 

where the index A labels the two independent transverse polarization states and according to asymptotic theory of 
interacting fields, Z is identified with the wave function renormalization for asymptotic states, namely Z is the wave- 
function rcnormalization constant in the vacuum. The constant Ck will be adjusted so as to subtract the photon 
number in the vacuum state. In free field theory Z = 1 , C/c = 1. In asymptotic theory the field creates a 

single particle state of momentum k with amplitude y/Z out of the vacuum state. The quantities Z , Ck account for 
renormalization aspects in the definition of the particle number in an interacting field theory. A detailed study in 
ref. |l5j shows that once Z and Ck are fixed by the wave function renormalization constant and the normal ordering 
subtraction in the vacuum, these terms cancel the zero temperature (vacuum) contributions. Therefore this definition 
of the number operator does indeed describe the number of photons in the medium at time t. 
Introducing the real-time correlation functions 

t')+g<{t,t')eit'-t)] (III.12) 
t)+g<{t,t')e{t~t')] (III.13) 

(III. 14) 
(III.15) 

where the superscripts ± refer to the forward (-I-) and backwards (— ) time branches. The expectation value of this 
number operator is defined in eq. (jll.43|l . and assuming A"^^ — Q;£^ t = expectation value in the time evolved 



{A+^^{k;t)A+,i-k:t'))=r^'\k) [g> it,t')e{t - 
{A-^^ik;t)A^,{-k;t'))^V=^\k) [g>it,t')eit' ~ 
{A-^^{k;t)A+^,{^k;t'))^r^^^{k) g>{t,t') 
{A^A~k-t')A+Ak-t)) ^V'^\k) g<{t,t') 
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density matrix is found to be given by 

dN 

x 

A 



cPk(Px 



(III.16) 



In terms of the center of mass field A introduced in eq. ()II.22|) it is straightforward to find that the correlation 
function in the bracket in eq. IjllLlljl is given by 



{Aimt^it')) = [g>{t,t')+g<{tX)] 

and the occupation number can be written in terms of the center of mass Wigner variable as follows 

1 



2kZ 



{A^,{t)-A_f:it)) + k'{Aj:it)-A_j:it)) 



-Ck 



(III.17) 



(III. 18) 



where the expectation values are obtained as in eq. (|II.43|I and A^{t) is the sohition of the Langcvin equation given 
by eq. (jllLlOp . 

Assuming for simplicity that = 0;£^ - = 0, the expectation value of the number operator eq. IjllLlljl as 



defined by eq. (|II.43p is given by 



^ A 



k Z \ 2k 

duj 
2^ 



ICTik,u;) [\J^kiio,t)\^ + k^\nk{u;,t)\' 



Ck 



(III.19) 



where we have introduced the auxiliary functions 



nk{LO,t) = / dTh{T)e-''^^ 
Jo 

Tk{Lo,t) = f drh{r)e-''^^ 
Jo 



(III.20) 
(III.21) 



and ICTikjU) is given by the fiuctuation-dissipation relation eq. (|11.58|l and Afk is the initial photon distribution 
function. 

This is one of the main results of this study. The expression piI.lQp is truly non-perturbative since the function 
fk{t) given by Ijlll.Sp describes the real time evolution after a Dyson resummation of the photon propagator in terms 
of the geometric series with the polarization S to lowest order in aem but in principle to all orders in a^. 

Thus it is clear that the formulation in terms of the non-equilibrium effective action for the physical transverse 
gauge fields described above leads to a novel non-perturbative approach to study photon production from a thermal 
source directly in real time. 

IV. LOWEST ORDER IN PERTURBATION THEORY: 

The results obtained above are general and as such may be unfamiliar within the context of photon production 
from a QGP in LTE. In order to establish the relationship to the usual approach we now obtain the photon yield 
in strict perturbation theory, namely to lowest order in aem- The Laplace transform (jIII.6|) to lowest order in the 
perturbative expansion in is given by 



fkis) = 



1 



T,T{k;s) 

s2 + P (s2 -|_ ^2)2 



(IV.l) 
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The fundamental solution can be readily obtained by inverting the Laplace transform using eq. (|II.60|) and is found 
to be 



(IV.2) 



with 



Sfkit) 



s'm[kt] 



— lmY.r^{k;uj)--^ -2-2 



t cos[A:t] 



did ImE^(/c; lu) 
TT 2k{u} — k) 



(IV.3) 



Inserting this perturbative solution in eq. piI.19|l and keeping terms consistently up to C'(Et) ~ 0{aem) and 
setting A/fc = we find to lowest order in aem 



(27r) 



3 dN{t) 
d^kd^x 









+ ^ r 




^ J-oc 



dw 



2ImS^(fc;w)n( 



1 — cos[(w — k)t] 

{OJ - fc)2 



(IV.4) 



The 1/Z in the first term is a consequence of introducing the wave function renormalization in the definition of 
the photon number in equation (|III.11|I . setting Z — l;Ck ~ I corresponding to the free field expression, the first 
term vanishes. The second term is exactly the one obtained in a perturbative expansion in real time in reference |l5j | . 
This equivalence can be inferred by noticing that the transverse part of the polarization IlT{k,Lu) as defined in eq. 
(IV. 9) in reference 15j and T,§:(k: w) as defined by eq. (|II.19p above imply that IlT{k,uj) = 2S^(fc;a;). Of course the 



perturbative study of reference [Igl defined the number of photons as in free field theory, corresponding to setting 
Z = l;Ck = 1. Thus we see that the perturbative evaluation of the expression (|III.18|I to lowest order in aem is 
exactly the same as obtained in the real time study in . 

Furthermore, as emphasized in ref.ll5l| the asymptotic long time limit takes the form 



1 



dN{t) 
d^xcPk ~ (27r)3 k 



2Imf:^(fc;w = k) 

eT — I 



+ 00 



duj 2 Iml].p(fc; u) 



1 



(IV.5) 



Therefore in the asymptotic long time limit and to lowest order in the photon production rate 



dN{t) _ 1 2lmf:§{k;iu = k) 



d^xd^k (27r)3 kieT-1) ■ 

is exactly the S-matrix result '15]. In the case that ImS^(fc;a; = k) = as for example in the hard thermal loop 
approximation|27n28|,i2^] a logarithmic dependence replaces the linear time growth in eq. ljIV.5|) 15.. 16.]. 

A. Photon production from kinetics: 

Kinetic theory provides an alternative approach to photon production. Within a simple kinetic description the time 
evolution of the number of photons in a phase space cell is given by a gain minus loss (master) type kinetic equation, 
which for the distribution function of each polarization is given by 

= [1 + n,,,(t)]r> - rH,A(i)r< (IV.7) 
with the photon number (assuming an isotropic distribution) per polarization 



(IV.8) 
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and ; are the forward and backward rates which are computed using Fermi's Golden rule which is equivalent 
to the S-matrix calculation of transition rates. Detailed balance entails 

rfc=e-^'rfc (IV.9) 

The solution of the kinetic equation H1V.7|I with the initial condition rife. a (0) ~ gives the following time evolution 
for the sum over the polarization 

nk{t) = 2n,g{k){l - e--^"*) ; 7fe = r< - r> (IV.IO) 

with neq{k) the equilibrium Bose-Einstein distribution function for photons. 

This solution illuminates two important aspects: i) the thermalization time scale is Tfc = and ii) for t « Tk 

nu{t) = 2n,q{k) 7fc t + 0{jI t") (IV.ll) 

the photon production rate is precisely determined by the expression for the photon number during time scales much 
shorter than that for thermalization, namely 

= 2r.,,(fc)7fe (IV. 12) 

The relaxation rate of the distribution function 7^. = 2rfe with Ffe relaxation (damping) rate of the single-particle 
excitation, in this case the photon damping rate, which is given by 

rfe = ^^^^?I%^ (IV.13) 

where the factor k in the denominator refers to the free photon mass shell. Accounting for the two transverse 
polarization state, one obtains the rate of photon production 

d Ukit) _ 2ImST(fc,cj = k) 



dt 



k{eT — 1) 



which is the result obtained above in eq. H1V.6(I . This simple kinetic description clearly shows the main physical 
assumptions and ingredients that enter in the computation of the photon production rate. 

The kinetic equation (jIV.7|) is obtained from a full quantum field theory description by defining the photon number 
operator as in the free field quantum theory of photons, thus the time evolution of this operator is solely determined 
by the interaction,24i|. The expectation value of the number operator is obtained in a perturbative expansion using 
the real time Feynman rules with free field photon propagators that include the distribution function just as in the 
equilibrium case'245. The main point of this discussion is that in the quantum field theory formulation leading to the 
simple kinetic equation (|IV.7|I the photons propagate as free particles, namely with the free photon dispersion relation. 
This ignores the fact that in the medium photons propagate as collective modes, not as single particle excitations. 
In order to include the collective effects in the medium, which are more relevant for photons with momenta smaller 
than or of the order of the temperature[23,|2l,|2^, a non-perturbative description is required. 

B. Kinetic interpretation of the S-matrix result: 

The analysis of the previous section leads to a kinetic interpretation of the S-matrix result. The S-matrix approach 
extracts the forward rate F^ by computing the transition probability in the infinite time limit, assuming n{k) = 
and ignoring the build up of the photon population, namely the terms with n^.x in the kinetic equation (|1V.7|I . thus 
leading to 



(IV.15) 
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Using the detailed balance relation eq. ljIV.9|) and the relation 7^ = 2Tk with Tk given by eq. (jIV.13|l one finds 
that eq. H1V.15|I is equivalent to the result H1V.14|I . 

Obviously the buildup of the population can only be neglected during a time scale i <C 1 /Tk , however the S- matrix 
approach manifestly takes the time to infinity extracting only the terms that grow linearly in time in this limit, 
and ignoring terms that grow slower and that can make important contributions during a finite time scale|l5l Il6j| . 
In assuming that the S-matrix approach is therefore valid only during a finite interval of time before the photon 
population builds up highlights the inconsistency of neglecting time dependent contributions associated with the 
finite lifetime. While this point has been emphasized in refs. 15, 16] the analysis based on the kinetic equation makes 
it explicit. 

A kinetic interpretation of the S-matrix result as gleaned from the full non-perturbative solution of the kinetic 
equation (|IV.10|I first assumes no initial population, that the rate is constant in time and that the population does 
not build up, leading to 



dNsM dNkin 



d'^xd^k d-^xd^k 
with the factor 2 accounting for the two polarization states. 



= 2r> 



(IV. 16) 



V. NON-PERTURBATIVE ASPECTS I: BREIT-WIGNER (NARROW WIDTH) APPROXIMATION: 



Having confirmed that the lowest order in the strict perturbative expansion in of the full expression ljIII.19p 
coincides with the results previously obtained in the literature, we now proceed to study non-perturbative aspects. 
The first stage of our study is to make contact with the kinetic approach to photon production studied in section 
HVH above. 

For weak couplings (aem o,nd as) when the width of the quasiparticle is much smaller than its energy, the photon 
spectral density 



p^(k]uj) 



lmt^{k-uj) 



A:2-RcS«(fc;a;) 



ImS^(fc, t 



(V.l) 



features a pole in the second (unphysical) Riemann sheet but near the real axis at the position of the "quasiparticle" 
pole, which is a solution of the equation 



w2(fc)-fc2-ReS^(fc;Wp(fc)) = (V.2) 

The imaginary part of the transverse photon self-energy evaluated a,i lj = ujp(k) determines the width or damping 
rate of the single quasiparticle excitation. If lmTjj^{k\u}p{k)) = then the quasiparticle is stable corresponding to a 
true pole in the physical sheet. If ImS^ (fc; ujp{k)) ^ the quasiparticle pole moves off the physical sheet becoming a 
resonance. In weak coupling the spectral density can be approximated by a Breit-Wigner form near the quasiparticle 
pole 

2 r 

p^{k; w w ujp{k)) = — (V.3) 

2nu)p[k) (w - Ljp(fc))^ + 

with the residue at the quasiparticle pole and the width given by 



= 1 



OReSg(fc;cj) 



ImEy (fc; u!p{k)) 
2ujp{k) 



(V.4) 



The residue at the quasiparticle pole Zk is in principle different from Z the wave function renormalization that 
enters in the asymptotic theory definition of the photon number (|lll.ll|l . The reason for the difference is that the 
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photons that are measured in the detector are asymptotic states, hence Z in this definition should refer to the vacuum 
value and not the in-medium residue at the quasiparticle pole. 

In strict perturbation theory the connection with the kinetic approach which applies to the distribution function 
of photons of dispersion relation ujp{k) = fc is established by taking Z — I] Zk — 1 ; Wp(fc) = fc ,C/t = 1, assuming that 
there arc no initial photons present, namely Afk — 0, and assuming the Breit-Wigner form of the spectral density 
(|V.3p in the full range of uj, not just near the quasiparticle pole. Under these assumptions 



1 Tfe _ Iml]^(fc;a; = k) 

2^ {uj - fc)2 + Tl ' " 2fc" 



P7ik;i^) = —— ,^,2 , T.2 ; Tfe- — (V.5) 



In the narrow width approximation ^ fc the fundamental solution is given by 



fk{t) 



sin[fci] 



(V.6) 



Inserting this solution and neglecting terms of 0{Tf,) 



in the numerator in eq. ljIII.19|) . and using the fluctuation- 



dissipation relation 



we find 



In the weak coupling limit, assuming a narrow Breit-Wigner spectral density, namely taking ImE|?(A:;a;) ~ 2krk 
and using that the spectral density is sharply peaked at lu ^ k we can replace = neq{k) with neq{k) the 

Bose-Einstein distribution function for photons. The integrals can be done straightforwardly and we find 



This result is the same as the solution of the kinetic equation ljIV.10|) since 7^ = 2Tk. The assumptions leading to eq. 
(|V.8|I . namely weak coupling, neglecting self-energy corrections to the dispersion relation as well as renormalization 
effects (wave-function renormalization) and the narrow-width-Breit-Wigner approximation for the spectral density 
are all approximations invoked in the Fermi's golden rule approach to the kinetic description. 

We note that even under all these assumptions, the result eq. (jV.81) is truly non-perturbative and a result of 
the non-perturbative Dyson- like resummation of the photon propagator manifest in the function fk{i)- The strict 
perturbative evaluation as presented in the previous section reveals that when Iml](fc,w ~ fc) 7^ 0, the photon yield 
grows linearly in time, a consequence of Fermi's golden rule (or secular terms). Such strict perturbative evaluation 
clearly is restricted to a time interval t <C FjT^ since for t 3> F^"'^ the photon yield attains the equilibrium form given 
by the Bose-Einstein distribution function. The secular terms (terms that grow in time) that appear in the strict 
perturbative expansion of eq. ljV.8|) are typically manifest as pinch singularities in finite temperature field theory [s^- 
A dynamical renormalization group program has been recently introduced that provides a resummation of these 
secular terms and leads to their exponentiation |3lj |. Thus, making contact with the dynamical renormalization group 
resummation of terms that grow in time, it is clear that the formulation presented above leads to a resummation of 
the perturbative series. The description of photon production given by our result eq. IjllFlOp not only coincides with 
the perturbative result in the strict perturbative expansion, but provides a systematic Dyson-like resummation of the 
perturbative series leading to a uniform expansion in the coupling valid at all times. This resummation is similar to 
that implied by the dynamical renormalization group approach of ref.[3l| and reveals the physics of the relaxation of 
the population towards equilibrium as well as the corresponding time scales. 



VI. NON-PERTURBATIVE ASPECTS II: PHOTONS OR PLASMONS? 



The discussion in the previous two sections above confirmed that the real time description of photon production 
describes both the lowest order results known previously and is capable of describing the thermalization of photons. 

We now study non-perturbative aspects of the production and propagation of photons in a locally equilibrated QGP. 

An important result from the hard thermal loop program in finite temperature gauge field theories|27l l28l [29| is 
that even for weak coupling the perturbative expansion must be resummed for momenta k < ujpi with ojpi the plasmon 
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mass. To lowest order in the HTL program for Nc = 3 and two (massless, up and down) quarks, the plasmon mass 
is given by 

/J 13.033 (^^) Mev (VI.l) 

which for example for the temperature expected at RHIC T ^ 300 Mev ujpi ^ 40 Mev. 

The real and imaginary parts of the transverse photon polarization in the HTL approximation are given by 
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(VI.3) 



In the HTL approximation, the position of the plasmon pole as a function of momentum k is determined by the 
solution(s) of the trascendental equation 
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The numerical solution for the dispersion relation is shown in fig. 
(transverse) plasma frequency ujpi — 39.099 Mev. 



(VI.4) 
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FIG. 1: The dispersion relation ujp{k) (in Mev) vs. k (in Mev/c) in the HTL approximation for uupi = 39.099 Mev, corresponding 
to a temperature T = 300 Mev . The solid line is the plasmon and the dashed line the free photon dispersion relations 
respectively. 



Even for momenta k 200 Mev ^ ujpi the difference between the plasmon and the free particle dispersion relation 
is about 15%. The real time solution fk{t) in the HTL approximation is given by 
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where is the residue at the plasmon pole (wavefunction renormalization) given by 
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The sum rule (|IIL9|I entails the following 
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(VI. 7) 



tj2-A:2-ReS«(fc;tj) + Inil]f(fc;w) 
The residue at the plasmon pole Zk is shown as a function of k for T ^ 300 Mev in figure |(5J) below. 



T=300 Mev 



500 1000 1500 2000 

k(Mev) 

FIG. 2: The residue at the plasmon (quasiparticle) pole for uipi = 39.099 Mev, corresponding to a temperature T = 300 Mev . 



Figure ^ below displays the continuum contribution fc{k,t) given by the integral expression in eq. (|VI.5|1 for 
the values k = 50; 100; 150 Mev/c respectively for T 300 Mev. It is clear from this figure that the continuum 
contribution is perturbatively small, of 0{aem) and damps out on time scales of 10 fm/c or longer even for large 
momenta. This time scale is of the order of the lifetime of the QGP expected at RHIC or LHC, hence fc{k,t) will 
contribute to 0{aem) during the lifetime of the QGP. 

The main point of this discussion is to highlight that once the photon is produced at a space-time vertex, it 
propagates in the medium as a plasmon collective mode, not as a free electromagnetic plane wave. The propagation 
of the photon in the medium as a collective mode has obviously nothing to do with the thermalization of the produced 
photon and only depends on the wavelength of the photon being smaller than the spatial size of the system. The 
situation is similar to an electromagnetic wave propagating in a dispersive medium with an index of refraction, the 
mean free path of the photon in the medium can be larger than the spatial extension of the medium itself, but the 
electromagnetic wave propagates (almost undamped) as a combination of normal modes in the medium. This point 
has been advanced by Weldon^^l within the context of dilepton production in a thermalized QGP^. 

This point can be simply illustrated by considering the average over the noise of the solution l|III.10l) with an initial 
value of the expectation value of the transverse gauge field and its time derivative corresponding to the case of a 



^ D.B. thanks Abhijit Majumder for discussions on this point 
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FIG. 3: The continuum contribution fc{k,t) given by the integral in JVI.5II for k = 150; 250; 500 Mev respectively as a function 
of t, for T ~ 300 Mev . 



positive energy plane wave, namely 



/2fc 



with J (A) being a transverse polarization vector. 

The propagation of this initial state is given by eq. ^TTT.injl . The noise average eq. (|II.38|I implies 



2V2k 
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Jp{k)t 



2V2fc 
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Wp(fc) 



.(VL9) 

Clearly the propagation of the transverse electromagnetic wave is not free. While the continuum contribution 
fc{k, t) vanishes at long time, its presence guarantees that the initial conditions (jVI.Sp are fulfilled via the sum rule 
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A. Beyond HTL 



The real and imaginary parts of the (transverse) photon self energy given by eqns. (|VI.2|1 . (|VI.3p are only valid in 
the hard thermal loop approximation k,u! <^T. The full one loop expression for the finite temperature contribution 
to the imaginary part of the transverse photon self energy is given by 



where ImE^^|,(fc;ti;) ; lm'S^2p{k;uj) are the Landau damping and two particle cut respectively, given bv[Tf 
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and the superscript (1) refers to one-loop. 

The real part of the self energy is found by the dispersion relation eq. (jll.51|l . The two particle cut contribution 
above the light cone would lead to an imaginary part for the plasmon since the plasmon dispersion relation lies above 
the light cone. This imaginary part is obviously related to the decay of the plasmon into hart massless fermion- 
antifermion pairs. Clearly this is unphysical, finite temperature self-energy corrections to the fermion lines lead to a 
thermal mass for the intermediate fermions which in the HTL approximation is at least 22% larger than the plasmon 
mass^, thus preventing plasmon decay into fermion quasiparticles. The correct imaginary part of the transverse 
photon self energy above the light cone requires a non-perturbative resummation that involves not only fermion self- 
energy corrections but also vertex corrections to satisfy the Ward identities. While the damping rate for a non-abelian 
plasmon at rest has been computed in ref. 33] and that for the abelian counterpart was estimated for large fc in ( 32]) 
the calculation of the imaginary part of the (transverse) photon self energy for all values of u, k is not yet available. 
What is available in the literature is the photon production rate calculated by the S-matrix approach to lowest order 
in aem and up to leading logarithmic order in ac|l2l]. From this result we can extract Im'Si^{u! — fc, fc) by making use 
of the relation eq. I|IV.6I) . 

The analysis presented above indicates that in order to obtain the complete real time dependence of fk{t), hence 
the real time photon yield, what is needed is the photon (transverse) self energy for all values of lo] fc. Below the light 
cone, the Landau damping contribution (jVI.ll|l is the leading order result, of 0{aemCtl) whereas above the light cone 
a resummation of the fermion lines and vertex is required for a consistent analysis of the plasmon width. Since the 
full imaginary part of the transverse photon self-energy above the light cone is not available, only its value on the free 
photon mass shell, namely lu = k, we model the imaginary part of the photon self energy in the full range of frequency 
as 



ImE(fc;w) w lmJ:^ll{k;uj) + lmf^{k;uj = fc) e{u}^ - k^) sign{uj) (VI.13) 

where ImE^^(fc; cu = k) is extracted from the S-matrix photon production rate eq. (jIV.6|l and the leading logarithmic 
result quoted in ref. [l^] , resulting in the following expression for the imaginary part on the free photon mass shell 



^ Even much larger than this estimate when the value for given by eq. IVl.lSi is included through gluon contributions to the fermion 
self energy. 
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The strong coupling as is a function of temperature, we will use the lattice parametrization |34| for the temperature 
dependence of the strong coupling (T) given by 
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(VI. 15) 



Although this lattice fit is valid at high temperatures and certainly not near the hadronization phase transition, we will 
assume its validity in the temperature range relevant for RHIC in order to obtain a numerical estimate of the photon 
yield. We note, however, that at the temperature expected at RHIC or LHC T 300 - 400 Mev as(0.3 Gev) 0.24 
and the validity of the perturbative expansion in as (even with leading logarithmic corrections) is at best questionable. 

The assumption that the imaginary part is constant above the light cone is consistent with the Breit-Wigner 
approximation. The one-loop Landau damping contribution below the light cone leads to the plasmon dispersion 
relation and is the leading order (0{aem)) contribution to the photon self-energy. The damping of the plasmon 
excitation is not correctly described by this Breit-Wigner form, since the correct plasmon damping rate is given by 
(see eq. HV.4p ) 
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However, since U0p{k) differs from k by less than 10% for the relevant range k > 200 Mev, and assuming that 
ImET(fc, Lu) is a smooth function of uj near tOpik) it is reasonable to assume that ImET(fc, uj — LUp(k)) ^ ImST(fc, w = fc). 
A reliable estimate of the error incurred with this approximation requires knowledge of IuiSt (k , co) for all values of 
UJ > k near the plasmon pole which is not yet available and its calculation is certainly beyond the scope of this article. 

Consistent with the above approximation for the imaginary part, the lowest order 0{aem) contribution for the real 
part is obtained b y th e dispersion relation eq. Ijll.51|l . In order to guarantee the vanishing of the magnetic mass in 
the abelian theory [23, we subtract the dispersion relation at zero frequency. The real part is therefore given by 
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where we have performed a subtraction at /cq = to ensure the vanishing of the magnetic mass and used the property 
(jll.53p . V stands for the principal part. 

The spectral density eq. (jV.ip with the real and imaginary parts of the self energy given by eq. I|VI.17|I and (jVI.lSp , 
(jVI.ll|l . (|VI.14p respectively features a sharp peak above the lightcone at a value of the frequency given by 
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We find numerically that the plasmon dispersion relation ujp{k) and residue at the plasmon pole Zk are remarkably 
similar to that obtained in the HTL approximation for momenta k > 200 Mev for 300 < T < 500 Mev displayed in 
figures H1I2|I respectively. This similarity was already pointed out in ref . |35| and our numerical study confirms the 
results obtained in that reference. This fact can be understood easily since for fc <C T the full one loop imaginary 
part reduces to that in the HTL approximation, and for k ^ T the corrections are truly perturbative. Above the 
light cone the photon spectral density can be very well approximated by the Breit-Wigner form 
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with the residue at the plasmon pole and the width given by 



lmT,^^{k;uj = k) 
2u)p(k) 



(VI.20) 



For 300 < T < 450 Mev ujp{k);Zk for k > 200 Mev computed from equations (|VI.13|) . HVI.17|) and (|VI.20ll are 
numericaUy indistinguishable from those shown in figures H1I2|I respectively, while is displayed in fig. Q). 
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FIG. 4: The approximation to the plasmon width given by 1VI.20I I for ujpi — 39.099 Mev, corresponding to a temperature 
T = 300 Mev . 



This figure in combination with the kinetic analysis of sections (|IV All reveal an important aspect: the solution of the 
kinetic equation (|IV.10|) shows that the time scale for thermalization of a photon with momentum k is tth;k = l/2rfc. 
Figure Q shows that for photons with k < 200 Mev the time scale for thermalization is < 15 fm/c which is of 
the same order as the lifetime of the QGP expected at RHIC and LHC. Hence low energy photons are expected to 
thermalize as quasiparticles in a QGP in LTE. This an important aspect that is completely missed by the S-matrix 
approach which would assign these photons a constant rate of emission even when they thermalize with the medium 
during the lifetime of the QGP. 

Below the light cone, the spectral density is determined by the one loop Landau damping contribution and is given 
by eq. IIV.l|l with the real part of the photon self-energy given by eq. ljVL17|) and the imaginary part given by 

t'il{k;Lu) given by eq. (IVl.lH . 

In this approximation, the real time solution fk (t) is given by 
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This solution represents a Dyson resummation of the lowest order contributions to the self-energy above and below 
the lightcone within the approximations detailed above. The initial condition fc{k, 0) = 1 is satisfied by the sum rule 
(|VI.7|I which we confirmed numerically. Thus the contribution from the Landau damping cut is necessary to fulfill 
the initial condition and the sum rule. 

With this real time solution, we now proceed to obtain the number of photons emitted from a QGP in LTE up to a 
time t by inserting this solution eq. (|VI.21|1 in the expressions l|III.19IIIL21|l . Since we are only computing the finite 
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temperature contributions to the photon self-energy and therefore neglecting the vacuum terms, we consistently set 
Z = 1 in eq. 1)111. 19fl . As discussed in detail in ref.|la| the zero temperature contribution to the self-energy yields the 
number of photons in the vacuum cloud which is precisely cancelled by the vacuum wave function renormalization. 
Consistent with this result in ref.fT^I we neglect the zero temperature contribution to the self-energy and set Z = 1, 
since both contributions cancel each other out in the photon number 15]. Thus we focus solely on the photons 
created in the medium. We have obtained the expression (|III.19ll under the assumption that the initial density 
matrix is factorized and gaussian for photons. Of course this assumption can be relaxed to contemplate initial states 
with correlations between quarks and photons, however as discussed in detail in ref.|l5j| such density matrix will 
not commute with Hqcd since any entanglement between quark and photon states will involve the quark current. 
Furthermore any correlated initial state will necessarily have to be modelled to include the initial preparation during 
and after the collision as discussed in detail in ref.|l5l|. In what follows we will study the simplest assumption, namely 
that the initial state is uncorrelated and with no photons. Accordingly, we set Afk = in eq. 1)111. 19|) . Introducing 
the solution l)V1.21|l into eq. ()III.19|I leads to many different contributions. We will explicitly keep terms that are of 
lowest order in aem in the uniform expansion in the electromagnetic coupling. We emphasize a uniform expansion 
to contrast with a naive expansion which features terms that grow in time, namely secular terms. In the many terms 
obtained we only keep those whose numerators are manifestly of 0{aem) and neglect those which are of higher orders. 
Specifically we find. 



= (VI.22) 
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with the wave function renormalization given to lowest order by 
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This expression summarizes the main result of this article: a non-perturbative expression the determines the direct 
photon yield directly in real time. It provides a Dyson-like resummation of the naive perturbative expansion, it is 
uniformly valid at all times and determines the yield to lowest order in agm- 

A strict perturbative expansion of this result up to 0{aem) leads to the following yield 
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Taking the asymptotic long time limit in this expression, the second term grows cx ln(t) at long times[TH hence it 
would be subleading and therefore neglected, while the first term yields the S-matrix result for the rate 
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with ImE^^ik^LO = k) given by eq. ()VI.14)) 
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Thus it is clear that our final result for the photon yield during a finite time t given by eqns. ljVI.22IVI.26p is 
a truly non-perturbative resummation of the naive perturbative expansion that incorporates real-time processes not 
captured by the S-matrix approach. 

Furthermore this result clearly reveals several physical aspects that are missed by the S-matrix approach and that 
underlie its shortcomings: 

• In the asymptotic long time limit, the only terms that survive are (|VI.23IVI.25IVI.26|I which in the perturbative 
limit Zk — I + 0{aem), reveal the correct physics: in a QGP in equilibrium with infinite lifetime, photons will 
thermalize as quasiparticles in the medium. Their abundance is determined by the Bose-Einstein distribution 
function for the correct quasiparticle frequency, and the wave function renormalization measures the weight of 
the quasiparticle in the spectral density. 

• The width Tk is formally of 0{aem), however expanding the expressions (|VI.23IVI.26p in Ffc, namely e~'"'=* ~ 
1 — Ffci + jF^t^ -|- 0{T^,t^) will lead to secular terms in time, namely terms that grow in time to all orders 
in perturbation theory. Similarly, si nce ujp {k) — k is formally of O(aem) one would be tempted to expand 
the argument of the cosine in eq. (lvT"2l . this again will lead to secular terms in time. Finally a naive 
perturbative expansion will attempt to replace LOp{k) — > k,Tk — *■ in the denominators of eq. ljVF24l) and 
(|VL26|) . however this will lead to logarithmic divergences at the Landau damping threshold which manifest 
themselves as logarithmic secular terms in time|l5j. The S-matrix approach computes the transition probability 
per unit space-time volume by taking the infinite time limit and extracts the linear term in time which gives 
a constant rate. Obviously taking the long time limit is manifestly ignoring the fact that however long the 
photon relaxation time, eventually photons will thermalize in the long time limit. A kinetic interpretation of 
the S-matrix result, as discussed in section (|IV Ap is given by 
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setting ujp{k) k, ^ 1 to lowest order in 0{aem), thus obtaining the result given by eq. (|IV.14p of section 



The non-perturbative aspects of the main result of this article, namely equations (|V1.22IVL26p make manifestl" 
clear the shortcomings of the S-matrix approach to extract the photon emission yield, including those in ref. 3 
such approach treats the plasma as a state in LTE of infinite lifetime and assumes that the photons are produced 
at a constant rate. These assumptions manifestly ignore the physically correct picture that if the lifetime of 
the QGP is infinite then photons produced will eventually thermalize. The validity of the S-matrix approach 
will then have to be justified on some intermediate time scale t <C 1/Ffc but then finite time effects must be 
included in the full calculation and the validity of the S-matrix approach will depend on the momentum of the 
photon. As we have discussed above, long-wavelength photons with k < 200 Mev thermalize on time scales 
of the order of the expected hfetime of the QGP at RHIC and LHC t ^ 10—20 fm/c (see fig. (0J), thus the 
S-matrix calculation is explicitly restricted to large momentum photons and only during a time scale t 1 /Tk 
which depends on the momentum. But clearly this is in manifest contradiction with the assumptions underlying 
the S-matrix calculation: the existence of asymptotic states and taking the infinite time limit. 

Assuming an initial state prepared in the asymptotic past not only ignores the true physical aspects of the 
problem, namely that quarks and gluons are partonic degrees of freedom in the incoming colliding hadrons, not 
asymptotic states, but also that however long the mean-free path of the photon, if the plasma lives long enough, 
these will thermalize. 



• In obtaining the final result given by equations ljVI.22IVI.26p we have systematically neglected terms that are 
truly perturbative and of higher order in an uniform expansion in aern- Namely the terms that we neglected are 
all bound in time and manifestly of higher order in aem- The final result (jVl.22IVl.26l) provides a resummation 
of the naive perturbative expansion much in the same manner as the dynamical renormalization group 31j which 
leads to correct kinetic eg nations (24]. 

Our main results and predictions for the yield and spectrum of the emitted photons are displayed below in figure 

© ,——.—1 

This figure shows the yield obtained from the result eqs. ljVI.22IVI.26p as compared to the S-matrix yield for 

T ~ 300 Mev;t = 10 fm/c which are values of the temperature and lifetime expected at RHIC. We have displayed 
the results only up to A: ~ 5 Gev since as discussed in |15J and above, for much larger momenta the wavelength of the 
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photon is likely probing scales shorter than the mean free path and the approximation of LTE is no longer reliable. 
Furthermore, as discussed in detail in [15iJ . for large momenta fc > 4 — 5 Gev the spectrum of photons is sensitive to 
the initial pre-equilibrium staged . 




FIG. 5: In (^4^) vs. fc(Mev) for T = 300 Mev and t = 10 fm/c. NRT{k) = dN/(fkcPx given by eqs. IIVI.22IVI.26ll and 
NSM(k) = dNsKi/cPkcPx given by eqs. HVI.29IVI.14fl . 



Clearly the real time yield is systematically larger than that predicted by the S-matrix approach. A flattening 
of the spectrum at fc ~ 2.5 — 3 Gev is a distinct feature that originates in processes that while subleading in the 
asymptotically long time limit dominate during the finite lifetime of the QGP. We have established numerically that 
amongst the four terms in eqs. I|VI.23IVI.26|) . the largest contributions to the yield during the QGP lifetime ~ 10 fm/c 
are those from eq. (|V1.23I) and cq. HV1.24|I while the contributions from ca. (|V1.25() . (|V1.26|1 are numerically smaller 
for the range of momenta displayed in fig. The first term given by l|VI.23(l features the fastest fall off with energy, 
while the remaining terms IIVI.24IVI.26p all contribute to the flattening of the spectrum. 

VII. DISCUSSION: 

A. Time evolution in URHIC: 

An important aspect of the real time formulation of photon production is the realization that particle production 
during ultrarelativistic heavy ion collisions requires the understanding of dynamics very different from usual collision 
experiments. In usual collision experiments the colliding particles in the initial beam are the 'in' states and the 
reaction products are the 'out' states. In and out states are eigenstates of the non-interacting Hamiltonian with the 
physical masses. The overlap between these states and those created by the Heisenberg field operators out of the exact 
vacuum at asymptotically early and late times is determined by the vacuum wavefunction renormalization constant 
Z. 

Time evolution in quantum mechanics and quantum field theory is an initial value problem and this is the basic 
ingredient in S-matrix theory approach to scattering. The S-matrix approach considers the time evolution of initial 
states corresponding to the 'in' states, namely eigenstates of the free Hamiltonian with the physical masses which are 
prepared in the infinite past. The S-matrix is the time evolution operator (in the interaction picture) from ti = —oo 
up to t/ = oo and evolves these initial states up to an asymptotically large time, at which the overlap between the 
'out' states, eigenstates of the free field Hamiltonian of physical mass is computed. The usual reduction formulas 
introduce the factor Z to account for the overlap between the states created by the full Heisenberg field operator and 
those of the free field creation operators of in and out fields. 

In URHIC the initial states are heavy ions and the final states are mesons and baryons, not the elementary quark 
and gluon fields in the QCD Hamiltonian. The collision involves the deconfinement, thermalization, and subsequent 
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re-hadronization of partons. The correct S-matrix approach should then include all aspects of these dynamical 
processes in order to provide a consistent framework for calculation. Needless to say, this is currently an impossible 
task. Instead, previous approaches to photon production seek to obtain the photon yield from the different stages 
independently. The pre-equilibrium stage is studied via parton cascade models which input parton scattering cross 
sections computed from S-matrix theory into a semiclassical transport formulation. The thcrmalizcd QGP state is 
studied by assuming a stationary QGP of infinite lifetime in thermal equilibrium with no memory of the initial pre- 
equilibrium stage. Photons from the hadronized state are studied by assuming a 'cocktail' of mesons. The results 
from S-matrix computations (which necessarily assume an infinite lifetime) in a thermalized QGP are combined with 
hydrodynamics to integrate the rate during the finite lifetime of the QGP, manifestly ignoring that the S-matrix 
approach has assumed an infinite lifetime. Photon production during the hadronized stage is calculated via S-matrix 
with in and out states corresponding to (short lived) mesons. Finally the total spectrum of photons is obtained by 
adding the contributions from the different stages. It should be obvious that there are glaring shortcomings in this 
program which justify deeper scrutiny. 

The real-time approach that we advocate is solidly based on the basic premise of quantum mechanics and quantum 
field theory, namely that time evolution is an initial value problem: given an initial state (or density matrix) and the 
Hamiltonian, time evolution is completely determined. Basic quantum mechanics states that given a density matrix 
at some initial time ti, p(ti), the density matrix at any time t (in the Schroedinger picture) is 

pit) = e-*^(*-*') p{t,) e'^(*-*') , (VII.l) 

where H is the complete Hamiltonian. In the case under consideration, H is the full Hamiltonian of QCD plus 
electromagnetism and is time independent. It is convenient to pass to the interaction picture of the non-interacting 
Hamiltonian Hq in which the evolution operator is given by the identity 

^-tH{t-U) ^ ^-iHo t jj^^^ ^^Ho U (VH.2) 

where U{t,ti) is the usual unitary time evolution operator in the interaction picture. 
The expectation value of any arbitrary Schrodinger picture operator is 

m» - (VIL3) 

Trp(t,) 

This is the same as the expectation value of the Heisenberg picture operator in the initial density matrix. Computing 
real time expectation values is very different from computing S-matrix elements. The former are in-in matrix elements 
as befits time evolution as an initial value problem, the latter gives a transition amplitude from a state far in the past 
to a state far in the future and is conceptually very different from obtaining an expectation value of an operator in a 
time evolved state. 

It remains to define the Heisenberg number operator: the Heisenberg field operator for photons creates a single 
photon state (in or out) out of the vacuum with probability ^/Z, with Z being the on-shell vacuum wave function 
renormalization constant. Since the number operator is a bilineal in the field it must be divided by Z in order to ensure 
that it counts asymptotic photon states with weight one, this is an important ingredient of the reduction formulae. 
Furthermore in free field theory a simple normal ordering suffices to ensure that the number operator annihilates the 
(bare) vacuum. However in an interacting theory, a normal ordering constant Ck must be introduced so that the 
expectation value of the Heisenberg field operator vanishes in the exact vacuum. This normal ordering constant can 
be simply extracted from the equal time limit of the operator product expansion (OPE) of the correlation function 
of the photon field. The short distance expansion is completely determined by the OPE in the vacuum. Therefore 
subtracting the constant Ck is necessary to define a normal ordered expectation value of the Heisenberg number 
operator consistently with the equal time limit of the OPE of the correlation function of the (transverse) gauge field. 

This is precisely the main point of our program, we obtain the expectation value of the photon number operator 
given by eq. (|III.11|1 which includes the overlap factor Z and the normal ordering subtraction Ck both determined in 
the vacuum consistently with asymptotic theory. In this article we implemented this program by obtaining the non- 
equilibrium effective action for the gauge field as explained in section l(ITj) . As was explicitly shown in our previous 
article (see |15)1. the vacuum wave function renormalization constant and the vacuum subtraction Ck cancel the 
vacuum contributions to the number operator. That this is indeed the number of photons (exact massless eigenstates 
of the free photon Hamiltonian) produced in the medium up to time t by the evolution is clear from the following 
trivial equality 

{Nk){t)^{Nkm+ f^^^^ (VII.4) 
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// the photons do not undergo further interactions after being produced during the Hfetime of the QGP, they will 
propagate freely towards the detector. 

As has been shown explicitly in section (|IV|) the S-matrix (or alternatively the kinetic approach) extracts the rate 
d < Nk > {t)/dt from Fermi's Golden rule, by taking the time interval to infinity which leads to a time independent 
rate, and inputs this rate in the calculation of the yield just as in eq. ljVII.4|) above accounting for hydrodynamical 
expansion. 

Of course if one knew how to evolve an initial heavy ion state through deconfinement, thermalization and finally re- 
hadronization with the full QCD Hamiltonian, one could begin with an initial density matrix of heavy ions at U —oo 
and compute the total number of photons produced all throughout the dynamics up to t — > +00 by implementing 
(|V11.3|I . Obviously this program is not feasible. However, the statement of formation of a thermalized QGP at a time 
of order ~ 1 fm/c suggests that at this initial time ti the density matrix is that of thermal equilibrium for the QGP. 
Time evolution being an initial value problem, if the initial density matrix is completely specified, then we can evolve 
it in time with the full time independent Hamiltonian H. Perturbation theory will be reliable if quarks and gluons are 
weakly interacting, the usual assumption for a QGT in LTE but certainly not during hadronization. Nevertheless, we 
can still obtain the number of photons produced up to any arbitrary time t from a weakly interacting QGP in LTE 
from eq. ljVH.3|) if we knew exactly what is the distribution of photons in the initial state. The main uncertainty in 
the calculation of the yield is precisely the knowledge of the photon distribution function in the initial density matrix. 

Only when photons reach thermal equilibrium with the medium will their spectrum and distribution be independent 
of the initial conditions. However if photons do not equilibrate, and this is the usual assumption, the spectrum of 
produced photons will depend on the initial state. In ref.j^ we have studied suitably modeled initial states and found 
that the photon spectrum is sensitive to the details of the initial state for energy larger than ~ 5 — 6Gev assuming a 
thermalization scale < Ifm/c. 

Recently a criticism has been leveled at the real time formulation for studying photon production from a transient 
QGplST^ that we propose. The authors in this reference claim that the unitary time evolution operator U{t, ti) in the 
interaction picture in eq. ljVH.2|l corresponds to 'switching on' the interaction at t — ti and 'switching-off' at time 
tf. This is clearly incorrect as the identity Eq. l|VH.2ll makes manifest. Such a switch-on-off of the interaction would 
entail a time dependent Hamiltonian with an artificial time dependent coupling. The identity (|VII.2p clearly states 
that time evolution is obtained with the full time independent Hamiltonian of QCD plus QED. This statement by the 
authors of ref. reflects a misunderstanding of the basic fact that time evolution in quantum mechanics and quantum 
field theory is an initial value problem: a state specified at some initial time ti is evolved in time up to an arbitrary 
time tf with the unitary time evolution operator given by eq. ljVH.2|) . As we mentioned above once the density 
matrix (or a pure state) is specified at any arbitrary time ti it can be evolved up to another arbitrary time t with the 
unitary time evolution operator. Any statements about 'switching on-off' of interactions is simply a misunderstanding 
of this basic fact in quantum mechanics. A series of formal steps are invoked in ref.|37j| to incorrectly argue that only 
taking tf 00 gives the total photon number. With this statement the authors are ignoring or glossing over a large 
body of work in the hydrodynamic approach to photon production from a thermalized QGP. In this approach an 
invariant rate, obtained assuming an infinite lifetime is input in the hydrodynamical evolution, and the total photon 
yield is obtained by integrating over the space-time history of the QGP, namely over a finite lifetime. One of our 
main criticisms of this approach is that the rate input in these calculations is obtained from an S-matrix computation 
which assumes an infinite lifetime, but is used during a finite lifetime. 

In a plasma that expands nearly at the speed of light, if its lifetime is infinite as the authors suggest, photons will 
ultimately thermalize with the plasma. Physically the QGP only 'lives' for a few fm/c, furthermore as the plasma 
expands and cools the coupling grows stronger and the theory becomes strongly coupled at the hadronization phase 
transition. Any formulation that attempts to use perturbation theory in a long time limit is obviously flawed. All of 
these physical arguments point out the limitations of an S-matrix approach and the inconsistency in taking an infinite 
time limit. 

The authors in ref. completely ignore the physical arguments that apply to the experimental situation therefore 
either ignoring or glossing over the inconsistency in taking an infinite time interval. Nor do they address the important 
issue of what is the initial state when the QGP is conjectured to thermalize, in particular what is the initial photon 
distribution function. S-matrix theory simply takes the initial time ti —00 and takes the initial states to be free 
field 'in' states. Clearly this is not the physical situation in a heavy ion collision and the pre-equilibrium stage is 
important. 

In ref. 37] the discussion stops short of actually calculating the number of photons in their limit tf 00. Had 
they done so they would have found a result that diverges linearly with t f . This is precisely the result of the S-matrix 
calculation where t/ is interpreted as the total 'reaction time' and the result is divided by to yield the rate. At this 
stage the physical question that the authors did not ask themselves is what is this 'reaction time'. Obviously if these 
photons are produced by a transient QGP, tf must be the lifetime of this state. Lest that the authors actually mean 
to follow the dynamics through the confinement and hadronization phase transition, in which case their perturbative 
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arguments are naive at best. 

The authors also question our definition of the Heisenberg number operator that includes the vacuum wave function 
renormalization Z and normal ordering constant C. As explained in the text and again explicitly in this section above, 
the wave function renormalization is required because the interpolating Heisenberg operator creates a single in or out 
photon state with amplitude \fZ out of the vacuum. This is a basic result of asymptotic theory and the factor XjZ 
is the usual factor that accompanies all bilinear correlation functions (such as the Green's function) in the reduction 
formulae, this is a basic issue in renormalization and in the LSZ formulation (see for example j38i|'). The normal 
ordering constant is a simple statement that the equal time limit in the operator product expansion gives a non- 
vanishing expectation value in the vacuum and must be subtracted. This is the basic statement that normal ordering 
and time evolution do not commute in an interacting field theory. Even in free field theory the number operator is 
defined with a normal ordering prescription (the usual symmetrization factor 1/2). Such normal ordering must be 
re-defined in the interacting theory and that is precisely what C achieves. It is defined at zero temperature, namely 
in the vacuum and as explicitly shown in ref.[l5j it cancels the vacuum contributions to the photon number in the 
interacting theory. 

In summary, the criticisms leveled in ref. have no credibility: i) statements such as 'switching-on' and 'off' reflect 
a misunderstanding of basic quantum mechanics and quantum field theory, ii) formal statements about infinite time 
limits etc, are only glossing over the physical situation under consideration, namely a transient state which undergoes 
a non- perturbative phase transition at a finite time, iii) the statements regarding wave function renormalization and 
subtractions also reflect a misunderstanding of basic renormalization issues in quantum fleld theory, iv) the authors do 
not even attempt to discuss the issue of the initial state, v) the criticisms in ref. [s^ do not attempt to scrutinize the 
usual assumptions with their glaring inconsistencies, nor do they provide a reliable framework to study the physical 
situation at hand, vi) the non-perturbative aspects associated with the propagation of photons in the medium (see 
discussion below) are not included in ref. [3^ , which our formulation in terms of the real time effective action describes 
consistently and systematically. Hence, these criticisms not only do not help to understand (or even address) the 
important and very relevant physical questions but disguise these fundamental aspects of URHIC physics. 



B. Non-perturbative aspects: 

The formulation of photon production based on the real-time effective action presented in this article allows to 
understand the following non-perturbative aspects of propagation and relaxation. 

i) Propagation: The photons produced in a QGP in LTE do not propagate as free waves in vacuum but in a 
medium with an 'index of refraction'. The propagation aspects of photons in this medium require that self-energy 
corrections to the propagator be included. A simple perturbative expansion of the propagator in terms of the photon 
self-energy E(a;, k) is 

G(a;,fc) = ^l-2+^^S(^,fc)^^ + --- (VII.5) 
LU^ — uj'^ — k'^ LU^ — k^ 

obviously features divergences on the photon mass shell, which in real time lead to terms that grow in time as for 
example featured in eq. ljVI.28|l . A Dyson (geometric) resummation of these self-energy corrections leads to the 
propagator 

^(-'^) = Z73p^ERfc) 

The poles in this resummed propagator describe the correct modes of propagation in the medium, namely the trans- 
verse plasmons and their relaxational dynamics. As usual the self-energy S(w, k) is calculated in perturbation theory 
up to a given order in the coupling, but the Dyson resummed propagator eq. ljVII.6p provides an all order resum- 
mation of select diagrams. It is precisely in this sense that the real-time effective action leads to a non-perturbative 
description which describes the correct propagation of photons in the medium as plasmon modes. In section IVll we 
discussed in detail the main aspects of photon propagation gleaned from the hard-thermal loop approximation to the 
self-energy. The poles in the Dyson resummed propagator determine the propagating modes, their frequencies and 
widths. Of course the dispersion relation and widths are obtained in perturbation theory up to the order in which 
the self-energy has been computed. This is the usual statement in many body or fleld theory that in order to obtain 
the position of the poles (propagating modes) a Dyson resummation of the self-energy corrections must be performed. 
These corrections that account for the propagation of photons in the medium with the correct dispersion relation are 
featured in the flnal equations (|VI.22|) - (|VI.26p in the wave function renormalization Zk, the dispersion relation ujp{k) 
and the width Tk in the denominators. 
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ii) Relaxation: Our final equation for the yield of photons produced during the lifetime of the QGP (jVI.22|l - 
(|V1.26|1 features contributions with the factor 1 — e"'"** with Fs. of 0{aem) and up to leading logarithms in a^. The 
naive perturbative expansion in aem replaces these factors by 1 — e"'"'"* ^ Ffct, which is the usual yield obtained 
from the S-matrix calculation. However, such approximation is obviously only valid for T^t <Si 1 (which again makes 
manifest the shortcoming of the S-matrix approach). However, as studied in detail in section FVI Al photons of energy 
~ 200Mev thermalize in the QGP during its lifetime of about lOfm/c and the product Tkt is a sizable fraction of 
unity for photons up to energies < 500 — 700Mev. Namely the full solution given by eas. (jVI.22|) - (|VI.26|) accounts 
for the partial thermalization of medium energy photons as a result of the non-perturbative Dyson resummation 
and is explicit in the exponentials which result from the complex poles in the Dyson resummed propagators. Naive 
perturbation theory certainly misses thermalization since the yield obtained from a perturbative S-matrix calculation 
will only feature a term of the form t and it also misses the correct propagation of photons in the medium. 

VIII. CONCLUSIONS, SUMMARY AND FURTHER QUESTIONS: 

In this article we have focused on studying production of photons from a QGP in LTE by implementing a non- 
perturbative formulation of the real-time evolution of an initial density matrix. The main ingredient is the real time 
effective action for the electromagnetic field exact to order and in principle to all orders in a^. The real time 
evolution is completely determined by the solution of a non-local stochastic Langevin equation which results in a 
Dyson-like resummation of the naive perturbative expansion of the photon self-energy. A quantum kinetic description 
emerges directly from this non-perturbative formulation. 

The main result for the direct photon yield from a QGP in LTE is given by eqns. (|III.19|I with (|III.8III.58p . 
We have confirmed that in a strict perturbative expansion, eg. I|III.19|I reproduces the S-matrix results as well as 
previously obtained results with initially uncorrelated states |l5l Il6t| . Furthermore we have shown that the result 
(|III.19|I reproduces the main features of quantum kinetics for the evolution of the photon distribution function. 

An explicit result for the photon self-energy to order and up to leading logarithmic order in ao|l2l| indicates 
that photons with momenta k < 200 Mev propagate as plasmon quasiparticles and thermalize on a time scale 
t < 10 — 15 fm/c, namely of the order of the lifetime of the QGP expected at RHIC and LHC. 

To leading order in a uniform expansion in aem and up to leading logarithmic order in as our main result is 
summarized by equations (|VI.22IVI.26|I and our prediction for the spectrum of photons from a QGP in LTE with a 
lifetime of order ~ 10 fm/c is displayed in fig. 

In the energy range of experimental relevance which we deem theoretically trustworthy 200 Mev ^ fc ^ 4 — 5 Gev 
the real time yield obtained in leading order and given by equations ljVI.22IVI.26|) is displayed in fig. |(SJ. It 
is systematically larger than previous estimates based on the S-matrix approach and features a flattening of the 
spectrum for fc > 2.5 — 3 Gev. We conclude that during the finite lifetime of a QGP in LTE expected at RHIC and 
LHC there are important real time processes that contribute to photon emission that are not captured by the usual 
approach. Furthermore, a distinct prediction resulting from our study is a flattening of the spectrum, a feature that 
is associated precisely with these real time processes. 

We have also discussed the limit of reliability of our results and the theoretical uncertainties that are present in any 
formulation of photon emission from a QGP in LTE: i) theoretical uncertainties in the initial state and pre-equilibrium 
stage, these will manifest themselves in the high energy part of the spectrum 13^ _15j, ii) a lack of knowledge of the 
imaginary part of the photon self-energy for all values of the frequency, iii) a likely breakdown of the reliability of the 
assumption of LTE for large energy photons. Once these uncertainties are resolved, a knowledge of the initial state 
(which itself requires understanding of the pre-equilibrium stage) and of the imaginary part of the photon self-energy 
can be systematically input in our non-perturbative formulation. 
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